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*0 Abstract 

i> 

We study the effects of an external magnetic field on the properties of the quasiparticle 
spectrum of the class of 2+1 dimensional strongly coupled theories holographically 
dual to charged AdS4 black holes at zero temperature. We uncover several interesting 
features. At certain values of the magnetic field, there are multiple quasiparticle peaks 
^ representing a novel level structure of the associated Fermi surfaces. Furthermore, 

increasing magnetic field deforms the dispersion characteristics of the quasiparticle 
peaks from non-Landau toward Landau behaviour. At a certain value of the magnetic 
field, just at the onset of Landau-like behaviour of the Fermi liquid, the quasiparticles 
and Fermi surface disappear. 



1 



1 Introduction 



With increasing interest in the various long-anticipated potential applications to condensed 
matter and nuclear physics of gauge/gravity duality techniques such as AdS/CFT and its 
generalisations [U El El Hj, it has been of considerable interest to understand how to cleanly 
capture the physics of the Fermi surface of strongly coupled systems using holographic meth- 
ods. This is because of the potential utility the methods might have to complement and 
extend the already impressive and ubiquitous power of the Fermi liquid paradigm, and also 
because there is a wealth of experimental phenomena that suggest that there are Fermi liq- 
uids that have behaviour that lies beyond the reach of standard techniques. Holographic 
methods may allow access to such behaviour. 

Recently, the work of ref. [5] presented convincing evidence for the Fermi surface of 
the strongly interacting 2+1 dimensional system at zero temperature and finite density that 
is holographically dual [51 E] to an extremal electrically charged Reissner-Nordstrom black 
hole in four dimensional anti-de Sitter spacetime (AdS^. (That work continued and refined 
the important initial studies presented in ref.jE], and there has been further discussion and 
extension of this particular line of development in refs. [21 HO] • See also ref. [TT] for a discussion 
of possible holographic descriptions of Fermi surfaces in a different context.) 

A key observation of ref. [5] was the fact that the Fermi liquid had distinctly non- 
Landau behaviour, showing unusual dispersion characteristics of the quasiparticle peak. This 
is perhaps to be expected, not just because the system is strongly coupled (remember that 
Landau Fermi-liquids arise in strongly coupled contexts too), but because one might an- 
ticipate that the effectively weakly coupled quasiparticle theory that arises (after strong 
coupling dressing) in the vicinity of the Fermi surface ought not to have a (simple) grav- 
itational dual. When all is said and done, it is, after all, just a free field theorj£| One 
might wonder, however, what range of dispersion characteristics might be accessible using 
holographic duals, and in particular how close one can get to Landau-like behaviour for a 
given system. The recent work of ref. JTU] answers some of this by doing a careful analysis of 
the AdS 2 x R 2 throat region region (first discussed in this holographic context in ref. [6]) that 
appears near the horizon of the black hole at zero temperature. Much of the physics of the 
critical exponents can be traced to the IR physics located down this throat and is controlled 
by the masses of fields there. 

Our work shows a different way of deforming the critical behaviour. As we will show, 

1 See, however, the work of ref. [9], but study it alongside the discussion in section VI. E. of ref. .10i. 



2 



adding a magnetic field moves the dispersion of the quasiparticle spectrum back toward 
Landau-like behaviour, in a manner that may have precise experimental analogues, especially 
given how natural magnetic fields are as laboratory probes and control parameters in a 
condensed matter context. Moreover, we find that at a given value of the magnetic field H 
(below what appears to be a certain limiting value H ma x), multiple quasiparticle peaks can 
appear in the spectrum, representing a finite series of levels (a kind of "band" structure). 
Our boundary condition is such that we are at the zeroth Landau level, so the levels we see 
are not the infinite family of equally spaced Landau-Rabi levels known from weak coupling 
intuition. 

To introduce a background magnetic field H we simply add magnetic charge to the 
AdS4-RN black hole, making it a dyon. We then study the system by probing it (as done in 
the H = case) with a minimally coupled Dirac fermion of electric charge q. We compute the 
retarded Green's function associated with this probe, from which, after Fourier transforming 
to momentum-frequency space (k,w), we extract our physics. 

Generically, it can be seen that the magnetic field shifts the effective mass of the 
fermion, which translates into a deformation of effective dimension of the operator we are 
probing the theory with. Following the discussion in refs.j9j ITU] , it is not hard to infer that 
this will certainly affect the dispersion of the resulting quasiparticle peaks, but this needs to 
be explored explicitly, and this is what we report on in this paper. 



2 Background 

The metric for the dyonic black hole in asymptotically AdS4 spacetime, using Cartesian 
coordinates is[T2"]: 

ds 2 = ^(-f(z)dt 2 + dx 2 + dy 2 ) + 1 ^^, (1) 
Z 1 v ' Z l J \z) 

F = 2Ha 2 dx A dy + 2Qadz A dt , 
f(z) = i+( H 2 + Q 2 )z i -(l + H 2 + Q 2 )z 3 = (l-z)(z 2 + z + l- (H 2 + Q 2 )z 3 ) . 

Here L is the length scale set by the negative cosmological constant A = —3 / L 2 in Einstein- 
Maxwell theory given by: 

^ = 2^ / d ' X ^{ R+ J2 - X F2 } ' (2) 
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where k\ = 87tGn is the gravitational coupling and our signature is ( — h ++)• The mass 
per unit volume and temperature of the black hole are: 

e = ^[l + Q 2 + H 2 } , T = ^[3-(Q 2 + tf 2 )] . (3) 

These are also the energy density and temperature of the 2+1 dual theory, which can be 
roughly thought of as living on the boundary at z — 0. The horizon of the hole is at z — 1. 
We choose our gauge such that the gauge field can be written in the following form: 

A t = 2Qa(z - 1) , A x = -2Ha 2 y . (4) 

The electric component sets, by virtue of its value on the boundary, a chemical potential 
\i = —2Qa for the U(l) charge, while the magnetic component determines a background 
magnetic field of magnitude T-L = —2Ha 2 .The parameter a has dimensions of inverse length. 
For simplicity, in the sequel we make the following redefinitions so that we are working 
entirely in terms of dimensionless fields and coordinates: 

t—>t/a, x^-x/a, y—ty/a, A t — >■ aA t . (5) 

2.1 The Probe Fermion 

We consider the Dirac action in this background: 

S D = J d A xyf^G i (*r M P M * - mM) , (6) 

where T>m is the covariant derivative given by: 

V M = d M + ^u abM T ab - iqA M , (7) 

and M are world indices and a, b are tangent-space indices. 0J a bM is the spin connection 
given by: 

(8) 

and 

T ab = 1 pfej ? T M = e M T a _ (g) 

The Dirac equation is given by: 

T M V M ^-m^ = 0. (10) 
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We choose \I/ such that: 

^ = 3/2.-1/4 -fcrf+k.* f (11) 

and 



(13) 



This reduces the equation of motion to: 

J— (d z =F m^/gZ) <f>± = ±z (7% + ry 2 ^ - 7 1 {2Hqy + fcj) 0=, 

V 9zz 



where: 



u=J^(u + 2 q Q(z-l)) , J g — = ^Tfi\ — = 4f ( 14 ) 

V -5« V 0** V -0« v/ 

Consider the limit of z — > of this equation. The solution is given by: 

+ = As" + (7 °" - y ^ ± K) + ^ B , (15) 

0„ = Bz~ m + ~* (7 °" + 71 (2ggZ/ + M - ^ A , (16) 

1 + 2m 

where A and B are functions independent of z. The various relationships between the 
various powers of z is relevant for the AdS / CFT dictionary [T3] . In particular, if we restrict 
ourselves to m > 0, then the source of the dual operator is proportional to B, and the vacuum 
expectation value (vev) of the operator is proportional to A. In this case, the dimension of 
the operator is given by: 

A = ^ + m. (17) 
We now proceed away from z = and further decompose the two-component fields <p± as: 

X± 

and make the following choices for the 7's: 

7 = ia 2 , 7 1 = ai , 7 2 = a 3 , (19) 
and these fields satisfy the following coupled equations of motion 



-(9 2 Tmv£)fe = T[-i< T + d y x T +i(2Hqy + k x )^ T ] , (20) 

Qzz 



-(4TmVfc)(± = T[iux T -d y £ T + i{2Hqy + k x )x T ] ■ (21) 

9zz 
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3 Specializations 



3.1 Zero Temperature 

Henceforth we will work at zero temperature. This means, from equation ^ that there is a 
relation between Q and H that we must bear in mind: 

Q = (3-H 2 y . (22) 

In particular, note that (we temporarily restore a for this discussion) the chemical potential 
\x = —2aQ naively gets shifted from its (T — 0, H — 0) value due to this relation. Also 
of note is the fact that as H increases, Q must decrease in order to keep the temperature 
vanishing. There is a maximum value of H, H m£LK = \/3, at which Q vanishes. This is all to 
be physically interpreted for our system as follows. First, we need to decide what physical 
quantity we are holding fixed while adding magnetic field to the system. A good such quantity 
is the chemical potential fi = —2aQ. Even though Q decreases for non-zero H, we simply 
increase a to hold it at the same value that it was at zero H. (According to equation (|3]), this 
increases the energy density e of the field theory, but this is quite natural as the magnetic 
field lifts the available energy levels of all charged particles.) Writing 3 — H 2 = Q 2 5 2 we can 
write a = a5~ 1 so that \i is fixed as 5 — > 0. Meanwhile we see that the physical magnetic 
field "H = —2a 2 H goes to infinity, showing that we do not have a limiting physical applied 
magnetic field. 

3.2 A Zero Temperature Ansatz 

In order to satisfy the appropriate boundary conditions at the event horizon, we redefine the 



fields as follows (for w^O; see subsection 3.4 for a comment on the u = case): 



X± = a ± (,,,)exp(^^J(l-^-^ 18 

^ A (A Ai(6<jrQ-4uO/18 



e ± = b ± (y,z)ex P ^^— ) j(l- Z y^-^^. (23) 

With this field redefinition, the equations of motion expanded at the event horizon give the 
following conditions: 

a+(y, 1) = b-{y, 1) , a.(y, 1) = -b+{y, 1) , (24) 
This result suggests that the "correct" variables to study the problem are actually: 

A+(y,z) = b-(y,z) - a+(y,z) , z) = -i (a_(y, z) + b + (y, z)) , (25) 
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B + (y,z) = a+(y,z) + b-(y,z) , £_(y, z) = i (b+(y, z) - o_(y, z)) . (26) 
The equations of motion are now given by: 
~Qxx ( r\ iw .-6qQ + 4w\ 

9 * + Tfx \? + N~ A + = ( 27 ) 



6(1 -z) 2 18(1 -z) 

fg^mB + - iuA + + i {d y B^ + (2Hqy + k x )B + ) , 



^r B .+ *±- +i -%j+ 4 » ]A _ = ,28) 



6(1 -z) 2 18(1-2) 

fg x ~ x mB_ - iuA_ - i (d y B + + (2Hqy + k x )B_) , 



6(1 -z) 2 18(1-2) 

- v fg^mA + + iuB + - i {d y A_ + {2Hqy + k x )A + ) , 

+ \? + z ~r^> r ) B - = ( 30 ) 



6(1 - z) 2 18(1 - 2) 

— y/g xx mA- + mB_ + z (d y A + + (2Hqx + fc^A-) . 

At the event horizon, we impose the following condition: 

d y B ± (y, 1) + 2HqyB T (x, 1) = , (31) 

which has solution given by: 

B+(y, 1) = B_(y, 1) = 5 exp (-Hqy 2 ) , 

B + (y,l) = -B4y,l)=B exp(Hqy 2 ) , 10 ) 

This solution can be understood as an infinite sum of the separable solutions of the problem 
in terms of the coordinate i] = y/2Hq (y + k x /2Hq). For this reason, we refer to this solution 
as the "infinite-sum" solution. We elaborate more on this relationship in our companion 
paper [T3j. This boundary condition restricts us to the zeroth Landau level. This then leaves 
us with the following conditions at the event horizon: 

A ± (r, 1) = , d z A±(r, 1) = t^j (K ± zm) B ± , (33) 

d z B±(x,l) = -^((48qQ-23uj)B ± Tl8VQ(k x Tim)d z A ± (x,l) 

Tl8V6d z (d y A T + 2HqyA ± )^ , 

i ( 54 \ 

= -W 48gQ-23w+ — (k 2 x + m 2 + AHqy (k x T tm)) ) B . (34) 
1U0 \ UJ J 
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Which solution is chosen depends on the direction of the magnetic field. For the rest of our 
analysis, we stick to the solution with B + = B_ which means we have negative H. 

In all that follows, we restrict ourselves to m = for our probe fermion. By the 
AdS/CFT dictionary, this means we are turning on operators with dimension A = 3/2. 
This case was studied in great detail for the case of H = in ref. [3j. 

3.3 Extracting the retarded Green's function 

Following the prescription for extracting the Green's function from ref. [13], the Green's 
function is given by (assuming m > 0): 

G * = !5S £ ( o J ' (35 > 

which in term of our fields A± and B± gives: 

r _ .X+(k) _ B + (k)-A + (k) .Uk) _. B-(k)-A.(k) 

+ ~ Uk)~ l B + {k) + A + {k) ~ \_(k) ~ 5_(Jfc) + A.(k) ■ W 

Note that these fields are given in momentum-space and not position-space. A subtlety 
here is that in ref. [13] , the fields always had a separable form with respect to the AdS radial 
coordinate and the space-time coordinates of the field theory, whereas here this is not the 
case. We assume that the prescription continues to hold and we simply Fourier transform 
our fields at the AdS boundary and use the result in the prescription. 

3.4 The case of to = 0. 

For the special case of u = 0, the condition A±(y, 1) = is no longer required. In fact, we 
need to change the behavior of the fields at the event horizon. Consider first the H = case. 



We go back to equations (23) and make the following redefinition: 



X± (z) = a ± (,)(l-,)^v^W ( (3?) 

t ± (z) = b±(z) (1 - ^V^W . (38) 

The exponent of the (1 — z) term is chosen such that there is a non-zero solution for the 
fields. Now we define A± and B± as before, and find the following condition at the event 
horizon: 

A+{1) = --=- (2qQ + »v^/3fc2 - 2qiQ*) B+(l) , (39) 
A-(l) = (2qQ + iV2,/3k 2 - 2g 2 Q 2 ) . (40) 
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In particular, if we take B + (l) = we have A_(l) = —A + (l), and we get for our 

Green's functions at the event horizon: 

Q ygg - 2q 2 Q 2 G _ ygg - 2g 2 Q 2 

+ V3k-V2qQ ' ~ v^ + v^gQ 
Another independent solution is to take: 

X± (*) = a ± (z)(l-z)-^ k2 ~^ Q2 , (42) 

= 6 ± (,)(l-,)--v^W (43) 

which gives at the event horizon: 

A+(l) = --j=- (2gQ - iV^fc 2 - 2g 2 Q 2 ) , (44) 

A_(l) = (2gQ - iy/2y/3k 2 - 2g 2 Q 2 ) . (45) 

Again if we take B+(l) = -B_(l), we have A_(l) = — A + (l) but with the Green's function 
gives as: 

q _ V3k 2 - 2q 2 Q 2 G ygg - 2g 2 Q 2 

+ V3k-V2~qQ ' V3k + V2qQ 

If we now turn on if, the situation becomes much more complicated. The equations of 
motion for A± and at the horizon require a more complicated exponent for (1 — z) in 
their behaviour. For our purposes, we will have no need to focus on the case of u = and 
determine these exponents since the interesting physics (specifically, the quasiparticle peaks) 
will appear away from ou — 0. We will therefore not pursue this further in the present work. 

4 Solutions 

4.1 A Separable Solution 

Recalling that we will work with m = henceforth, we begin by considering the trivial case 
of a separable solution which occurs for k x = 0, i.e. 

A±(y, z) = A(z) exp (-Hqy 2 ) , B ± (y, z) = B(z) exp (-Hqy 2 ) . (47) 

In this case, the four equations of motion reduce to two: 



/ 9xx 


(d z A- 


iuj 


9zz 


6(1 -zf 


9 XX / 




iu 


9zz 


6(1 -z) 2 



18(1 -z) 
-6qQ + Au 
18(1 -z) 



9 



Therefore in the separable case, Gr(1, 1) = Gr(2, 2). In fact, looking at the initial conditions 
at the event horizon for A(z), we find that d z A(l) = A(l) = 0, and therefore we can 
consistently take A(z) = 1. Therefore, we automatically get that Gr(1, 1) = i. 
The other separable solution occurs for H = 0, and our solution should match those of 
ref. [SJ. In this case, we can take: 



At (y,z) = At (z) , B ± (y,z) = B± (z) . 



(50) 



Choosing q = —0.5 we do indeed recover the results of ref. [SJ, as we show in figure [TJ 
These two cases of separable solutions provide good test cases for the accuracy of numerical 



Im(G R (2,2)) 




Figure 1: Solution for H = and k x = 1.2 using Mathematical NDSolve 



procedure when we solve the infinite-sum cases next. 



Before proceeding, we emphasize that our separable ansatz in equation (47) only 



considers the zeroth order Hermite function. The ansatz can be generalized for higher order 
Hermite functions, the physics of which we study in a companion paper |14j . 



4.2 Infinite— sum cases: Full Numerical Method 

We now proceed to solve the problem with k ^ and H ^ 0. We choose to discretize the 
equations of motion using a "fully-implicit" method: 



A + (yj,z n ) = A + (y j ,z n -. 1 )-iAzF + (z nr .i)A + (yj,z n - 1 ] 



(51) 



+i- 



Az f B_(y j+1 



\ffiZn-l) 

where we have defined: 



B_(y j _ 1 ,z n _ 1 ] 



2Ay 



+ {2Hqyj + k x ) B + (yj, z n . 



Fjz) 



CO 



-6qQ + Au 



u 



Q(l-zf 18(1-*) vT^J 



(52) 
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The problem is now to solve a large number of linear equations at each z step. The associated 
matrix is of block tri-diagonal form, i.e. : 



\ 



A 2 




Ci 
A, 





c 2 

B 3 






c 3 



V 



(53) 



/ 



where each of the Ai, B iy Cj's are matrices (8x8 matrices in this case). This matrix can be 
inverted using a "forward elimination-backward substitution" method. 

4.2.1 Accuracy of Numerics: Two Tests 

The first test of our numerical method is to reproduce the known separable results we describe 
above. We present the result for H = 0, (with k y = and non-zero k x ) in figures [2] and [3j 
demonstrating with the latter that we can reproduce the correct quasiparticle peak and its 
associate scaling behaviour. It is interesting to test the case of k x = 0. 



Im(G fi (2,2)) 

7 r 




\m(G R (2,2)) 

It 




(b) k x = 3 

Figure 2: Solution for H = using implicit method to solve PDE. 



We find considerable deviation from the expected answer of Gr(2, 2) = i. We can 
understand this completely as arising from numerical error. The field A + at the boundary 
should be zero in this case, but it turns out to be proportional to the derivative of the B's. 
It is straight-forward to understand this discrepancy: it comes from the finite difference 
approximation to the ^-derivative, which has zeros for non-zero H at zero k x , causing a 
loss of accuracy. Given that we have rotational symmetry, we can choose, without loss of 
generality, to put all our momentum into k x and work in the more well-behaved k y = 
sector for the remainder of our investigations. 
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We are now ready to numerically explore (k, u) landscape for various H, searching 
for quasiparticle poles in the Green's function. 



5 Observations 

There are several important general features that we encountered in our numerical explo- 
rations, and we collect them all together here. 



5.1 A Gap, A Ridge, and Some Poles 

• For a given magnetic field, as k x is increased, a region appears where Im(G#(2,2)) is 
negative for values of uj < u; gap . We explain the meaning of this region below, but 
we illustrate in figure [4] how u gap behaves as we increase k x for fixed magnetic field. 
In particular, we note that as k x — > oo, the behavior of u gauP is linear and appears to 
remain a fixed distance from the u = k x line. For higher magnetic fields, the separation 
from the line decreases. 

• We find a peak for values of lo > u gap . We refer to the position of the peak as a;*. We 
show samples of the behavior of to* in figure [5] This is the ridge of ref. [J)J. From figure 
|5j one sees that as the magnetic field is increased, the ridge gets closer to the w = k x 



line, crosses it, and continues to approach it from above (see figure |5(b) ) . In addition 



the ridge always remains above the gap line (except when they meet). We show this 
in figure [6j We also show an example of the large k x behavior for different magnetic 
fields in figure [7| 
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• For a fixed value of H, when u gap = at some k x = kp, we find a pole in ImG#(2, 2). 
Notice that this means that a pole appears whenever the ridge curve touches the gap 
curve. We find that the curves never cross but bounce off each other. We show 
an example of such a pole in figure |8j This is a quasiparticle peak. It indicates a 
Fermi surface at kp. Away from the poles, the numerical results should be interpreted 
carefully since there one should consider a complex u and not a purely real u. This is 
relevant because in our numerics the poles occur immediately after a region where the 
imaginary part of the Green function is negative, which is not allowed if the theory is 
unitary. We believe that the proper treatment of the complex u in these regions would 
resolve this. 

The fact that u > tells us that the Fermi energy Ep is greater than the U(l) chemical 
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0.05 0.10 0.15 0.20 0.25 0.30 0.35 



Fi gure 7: Behavior of w ga p (red) and (blue) at fixed k x = 100 

potential //. This is quite natural: the physical magnetic field H = —2a 2 H has lifted 
the energy of the charged particle, as earlier discussed. (Note that we see peaks for 
negative u as well, but focus on positive u for physical interpretation.) 

In particular, we note that the behavior is very reminiscent of the behavior of the full 
Green's function (that includes both the advanced and retarded Green's function) at 
the Fermi surface where both quasiparticles and quasiholes are observed. It is also 
important to note that there can be multiple (but we believe only a finite number) 
such peaks, at a given H, appearing at distinct values (u*, k x ). 

• There exists a maximum magnitude of the magnetic field at approximately \H\ > 0.37 
above which we do not find any poles. We notice that above this maximum, the ridge 
remains above the w = k x line for all values of k x and the gap remains below the line, 
preventing the two from meeting. We explain in the following subsection what this 
maximum corresponds to. 
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Im(G 22 ) 




Figure 8: Behavior of the pole in the TmG R (2, 2) for H = -0.1 at k x =0.871847 

5.2 Dispersion of the Quasiparticle Peaks 

We proceed to study the scaling properties of the poles observed in Im(G fi (2, 2). We present 
samples of the results in figure [9] We find that the behavior near the pole seems very close to 
linear once H ^ 0, if approaching from below or above. Contrast this with the case of H = 
(discovered in ref. [5], recomputed here and displayed in figure [3J. However, as one can see 
from the figure, there is always a change in slope as we go from k < kp to k > kp. This 
suggests that there is perhaps curvature to the lines infinitesimally close to kp. Therefore, 
the system is still in a non-Landau Fermi liquid regime but with a scaling less than that 
found for the H = case. 

Crucially, we see that as we increase the magnetic field, the slopes of the two lines 
begin to approach each other. Our results suggest that when the slopes approach the limit 
of being equal, that is exactly when we no longer find any poles. This seems to be the 
mechanism by which the system protects itself from becoming a simple (linear) Landau 
Fermi liquid, as discussed in the introduction. 

In addition, we note that when there are are multiple poles present (for a given H), 
the dispersion behavior of the poles may not be the same. We present an example of this 
in figure [TUj The results suggest that the further along the ridge one finds a pole, the more 
the behavior of the pole approaches that of a Landau Fermi liquid (i.e., linear dispersion). 
This suggests why there may only be a finite number of poles for a fixed magnetic field, since 
as we saw earlier, beyond a certain H, or far enough along the ridge, the ridge and gap no 
longer cross. 

6 Conclusions 

We've found a rich set of physical results (the deformation of dispersion characteristics, and 
the discrete Fermi levels) from our holographic studies of the Fermi surface and quasiparticle 
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Figure 9: Scaling behavior of W*. The red dot is where the pole is located, which corresponds to 
where the two lines intersect. 



spectrum in a background magnetic field. We expect that there is even more rich physics to 
be found from these systems. It is very exciting that some of these phenomena seem akin 
to the sorts of strongly coupled physics that are experimentally accessible, including with 
background magnetic fields ClS db probe of physics. 

Our exploration work was primarily numerical. A precise analytic relation between 
the value of the magnetic field and the nature of the dispersion of the peaks would be a 
valuable result. However, there is an apparent obstruction to doing the obvious generalization 
of the analysis of ref. [TU] that exploits the presence of the AdS2 x IR 2 region. There, an 
expansion about the uj = point, where the quasiparticle peak appears, representing a 
Fermi surface with Fermi energy Ep equal to the U(l) chemical potential fi. 

Things are different in our case. While it is trivial to show that there is again an 
AdS2 x IR 2 region near the horizon of the dyon (rather nicely, the electric field is entirely 
in the AdS2 and the magnetic field is entirely threading the IR 2 ), this is not enough. In the 
presence of magnetic field, the energy of the system gets lifted. The lowest available level 
at which we find a quasiparticle peak, for a given T-L has Fermi energy Ep greater than the 
£7(1) chemical potential /i, and so occurs away from the lo = point that was the focus of 
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Figure 10: Scaling behavior of for two poles at fixed magnetic field H = —0.35. The red dot is 
where the pole is located, which corresponds to the two lines intersect. 



refs.[5l HH]- There is no good analytical guide (so far) as to where in u the peaks will arise, 
and so we search for them numerically. An analytical characterization of exactly how the 
presence of % affects the dispersion relation will have to await future work. 
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